Experimental Realization of a 2D Fractional Quantum Spin Liquid 
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The ground-state ordering and dynamics of the two-dimensional (2D) 5=1/2 frustrated Heisenberg 
antiferromagnet CS2CUCI4 is explored using neutron scattering in high magnetic fields. We find 
that the dynamic correlations show a highly dispersive continuum of excited states, characteristic 
of the RVB state, arising from pairs of 5=1/2 spinons. Quantum renormalization factors for the 
excitation energies (1.65) and incommensuration (0.56) are large. 



The concept of fractional quantum states is central to 
the modern theory of strongly correlated systems. In 
magnetism, the most famous example is the spin S=l/2 
ID Heisenberg antiferromagnetic chain (HAFC) where 
pairs of 5=1/2 spinons are deconfined from locally al- 
lowed S=l states; a phenomenon that is now well es- 
tablished both theoretically |jl|] and experimentally 
These spinons are topological excitations identified with 
quantum domain walls. Experimentally, such fractional- 
ization is manifest as a highly dispersive continuum in 
the dynamical magnetic susceptibility measured by e.g. 
neutron scattering ||^, and for the HAFC identified as 
creation of pairs of spinons. 

In 1973 Anderson ^ suggested that a 2D fractional 
quantum spin liquid may take the form of a "resonating 
valence bond" (RVB) state comprising singlet spin pair- 
ings in the ground state, and with pairs of excited 5=1/2 
spinons separating via rearrangement of those bonds. 
The dominant feature of the RVB state, present in all 
its theoretical descriptions [^-|6j is an extended, highly- 
dispersive, continuum. To date this feature remains un- 
observed in any 2D magnet; in the case of the 5=1/2 
Heisenberg square lattice (HSL) mean field confining ef- 
fects lead to 5=1 magnons and a renormalized classical 
picture of fiuctuations around local Neel order emerges 
One may think, however, that because frustrating 
interactions can counteract the staggered fields respon- 
sible for confinement |^,||, they may provide a route to 
generating fractional phases in 2D. 

We explore such a scenario by making neutron scatter- 
ing studies on CS2CUCI4. By exploiting its unique exper- 
imental properties as a low-exchange quantum magnet 
we reveal an unexpectedly strong two-dimensionality 
in the form of a triangular antiferromagnet with partially 
released frustration. The simplicity of the couplings in 
CS2CUCI4 makes it a model system to investigate generic 
features of 2D frustrated quantum antiferromagnets. 

The structure of CS2CUCI4 is orthorhombic {Pnma) 
with lattice parameters a=9.65 A, &=7.48 A and c= 12.35 



A at 0.3 K. Magnetic interactions are mostly restricted 
between Cu^+ 5=1/2 spin-sites in the {b,c) plane, see 
Fig. 0(a), with coupling J along b ("chains") and zig- 
zag "interchain" coupling J' along the c-axis A 
small interlayer coupling J" < 10~^ J (along a) stabilizes 
3D order below T/v = 0.62 K into an incommensurate 
structure along b due to the frustrated couplings; weak 
anisotropics confine the ordered moments to rotate in 
cycloids near coincident with the (6, c) plane, see Fig. 
|l|(d) but with a small tilt of the cycloidal plane relative 
to the (6, c) plane whose sense alternates along c such 
that for each plaquette (isosceles triangle) (Si_(S2 x S3)) 
is small but nonzero (order is noncoplanar) |ll| , making 
this system a candidate for a chiral spin state jlj]. The 
minimal Hamiltonian determining the magnetic order is 

with each interacting spin-pair counted once, see Fig. 
|l|(a); a detailed description of the full Hamiltonian in- 
cluding small Dzyaloshinskii-Moriya terms is given else- 
where |l3|] . The Hamiltonian interpolates between non- 
interacting HAFCs (J' = 0), the fully frustrated trian- 
gular lattice (J' = J), and unfrustrated HSL (J = 0). 

Quantifying the couplings in (1) is of considerable im- 
portance both to guide theory and put our results in con- 
text. We do this using the following approach: neutron 
diffraction measurements were made on a single crystal of 
CS2CUCI4 in magnetic fields up to 7 T and temperatures 
down to 0.2 K using the PRISMA time-of-flight (TOF) 
spectrometer at the ISIS spallation neutron source. For 
fields along a (near perpendicular to the planes of spin 
rotation) a 3D incommensurate "cone" order is stable 
up to full ferromagnetic (F) alignment {Be — 8.44 T at 
T=0.03 K), see Fig. |l|(c). At T=0.2 K magnetic Bragg 
peaks arising from the transverse spin rotation move from 
£0=0.030(2) in zero field to e=0.047(2) at 7 T, where e 
is the incommensuration relative to Neel order, see Fig. 
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|l|(f). Mean-field theory predicts no change with field, and 
the large renormalization observed is a purely quantum 
effect. Since the ferromagnetic (F) state is an eigenstate 
of (|l|) with no fluctuations, e at the saturation field Be 
is at its classical value [Q siuTrec = J'/2J. Higher- 
field measurements Q observe ec = 0.053(1), imply- 
ing an exchange coupling ratio of J'/ J=0.33(l). The 
resulting quantum renormalization of the zero-field in- 
commensuration eo/ec=0.56(2) is similar to the predicted 
value of 0.43(1) (J'/ J=0.33(l)) estimated by series ex- 
pansions using a paired singlet basis [ p^ Additionally, 
the determined exchange coupling ratio is in agreement 
with the observed 2D dispersion in the saturated phase 
at 12 T II a which give the bare exchange couphngs- 
per-site J=0. 375(5) meV within chains and 2J'=0.25(1) 
meV between chains. This demonstrates that "inter- 
chain" couplings are of the same order as "intrachain" 
and Cs2 CuCli is therefore a quasi-2D system. These ob- 
servations require a change in the point of view taken 
by earlier studies [0, which proposed a quasi- ID pic- 
ture based on estimates of the "interchain" couplings not 
including the large quantum renormalization of the in- 
commensuration reported here. We now present detailed 
measurements of the excited states. 

Dynamical correlations in a 2.5 cm'^ single crystal 
of CS2CUCI4 were probed using the indirect-geometry 
TOF spectrometer IRIS, also at ISIS. The energy resolu- 
tion (0.016 meV full-width-half-maximum (FWHM)) was 
an improvement of nearly an order of magnitude com- 
pared to our previous studies [Q. The detectors form 
a semicircular 51-element array covering a wide range 
of scattering angles (25.75° to 158.0°). The sample was 
mounted with the (a, b) scattering plane horizontal in a 
dilution refrigerator insert with base temperature 0.1 K. 

The results for the b* dispersion are shown in Fig. 
§(a). Although the data points also have a finite 
wavevector component along a*, no measurable disper- 
sion could be detected along this direction confirming 
that the coupling between layers J" is negligible. The 
observed dispersion is well accounted for by the prin- 
ci pal spin-wave mode |l5| of Hamilto nian (0), w(k) — 
\/ (-'"k - -'"q) [(-'"k-Q + Jk+q) /2 - Jq] where the Fourier 
transform of exchange couplings is Jk = J cos 27rfc -|- 
2 J' cos tt/c cos t:1 and k = (h, fc, I). The ordering wavevec- 
tor in the 2D Brillouin zone of the triangular lattice is 
Q=(0.5 -I- eo)b* and the effective exchange parameters 
J=0.62(l) meV and J' = 2Jsin7reo (fixed) are in agree- 
ment with |ic| ] . The quantum renormalization of the ex- 
citation energy J/J=1.65 is very large and is similar to 
the exact result 7r/2 for the ID 5=1/2 HAFC (see e.g. 
0). In contrast, the spin- wave velocity (energy) renor- 
malization in the unfrustrated S'=l/2 HSL is only 1.18. 
Such large renormalizations of energy (1.65) and incom- 
mensuration (0.56) show the crucial importance of quan- 
tum fluctuations in the low- field state of CS2CUCI4. 



A remarkable feature of the measured dynamical corre- 
lations is that these do not show single particle poles, but 
rather extended continua. Fig. ||(b) (open circles) shows 
a scan at the magnetic zone boundary taken at 0.1 K. 
The scattering is highly asymmetric with a significant 
high-energy tail. The non-magnetic background (dashed 
line) is modeled by a constant-plus-exponential function. 
The magnetic peak disappears at 15 K (solid circles) and 
is replaced by a broad, overdamped, paramagnetic sig- 
nal. Figs. I^A-D show 0.1 K data properly normalized 
and corrected for absorption, and with the non-magnetic 
background subtracted. To quantify discussion of the 
dynamical correlations measured by neutron scattering, 
we first consider a spin-wave model, which is known to 
provide a good description of the unfrustrated HSL @] . 

The dynamical correlations of the spin-wave model 
|]l5| for Hamiltonian (|l|) exhibit single-particle poles from 
three spin-1 magnon modes, polarized with respect to 
the cycloidal plane. Fig. ||(a) shows the main dispersion 
mode i-u{k), polarized out-of-plane, and the two secondary 
modes, uj^{k) — uj{k — Q) and uj^{k) = uj{k + Q), both 
polarized in-plane, where the equilibrium spin direction 
rotates in-plane with wavevector Q. The expected scat- 
tering is given by the dashed lines in Figs. ^A-D, which 
clearly fails to account for the observed intensity as well 
as the extended high-energy tail of the scattering. This 
tail is not an instrumental effect as the FWHM of the en- 
ergy resolution (horizontal bar in Fig. |^(b)) is an order 
of magnitude narrower than the signal width. Including 
next-order processes also fails to account for the scat- 
tering: the two-magnon scattering (polarized in-plane) 
contribution to the lineshape is also shown in Fig. |3|D 
(shaded area) - this was calculated numerically using the 
method described in | p6[ | using the experimentally esti- 
mated spin reduction AS* ^ 0.13 [0 to normalize the 
elastic, one- and two-magnon scattering. 

Because in a neutron scattering process the total 
spin changes by AStotai—0,a the absence of single- 
particle poles and the presence of excitation con- 
tinua implies that the underlying excitations carry frac- 
tional quantum numbers. For J'—O, these are rig- 
orously known to be 5=1/2 spinous 0, and two- 
spinon production is the principal neutron scattering 
process |2j. Our analysis shows that the measured 
scattering can be described by the Miiller ansatz line- 
shape appropriate to the ID J' = limit S{k,uj) ^ 
8 (w — LOi{k)) 6 (wtt(fc) — w) I \/l>j^ — oj'i{k), where Q is 
the Heaviside step function, and uji and cou are the lower 
and upper continuum boundaries, generalized to 2D such 
that: (1) the total cross-section has three continua with 
the lower boundaries uj{k), Lo~(k) and a>+(fc) shown in 
Fig. |^(a); (2) the continua have equal weights and are 
isotropic in spin space; and (3) a modified upper bound- 
ary LOu (dashed upper line in Fig. ^(a)) is used. This 
model provides an excellent description of the data, see 
Figs. P|A-D. It is also noteworthy that both the asymmet- 
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ric dispersion, characteristic of 2D frustrated couplings, 
and the excitation continua are essentially unchanged at 
T—0.9 K above r;v=0.62 K in the disordered spin liquid 
phase showing that ordering affects only the low-energy 
behavior. We conclude that, in contrast to the HSL 
(J=0) where unfrustrated couplings confine spinons into 
S~l magnons throughout the Brillouin zone, CS2CUCI4 
has fractional spin quasiparticles carrying the same quan- 
tum numbers as in the HAFC (J'^0), namely S=l/2 
spinons, and, further, that these spinons are modified by 
the two- dimensionality at all energy scales. Although no 
evidence of spinon confinement is observed at any of the 
energy scales probed in our experiments, low energy S=l 
Goldstone modes are expected to occur in the 3D ordered 
phase; weakly coupled HAFCs have recently been shown 
to exhibit dimensional crossover in the dynamical cor- 
relations from low-energy 3D spinwaves to high-energy 
ID spinon continua on an energy scale of the interchain 
couphng [0. 

Magnetic fields applied within the ordering plane have 
a profoundly different effect from those along a. In fields 
along c a transition occurs above 1.4 T (T <0.3 K) to a 
phase, marked S on Fig. 1(b), where (1) the structure 
is elliptical with a large elongation along the field direc- 
tion, (2) the incommensuration approaches a linear rela- 
tion with field with a large slope [|l^, see Fig. 1(e), and 
(3) the total ordered moment decreases with increasing 
field. Above 2.1 T there is no long-range order at least 
down to 35 mK, and the system is in a spin liquid state. 
In this phase the dynamical correlations show shifts of 
continua and redistribution of scattering weight pO| ] as 
expected for spinon states and compared in pO[ | with ID 
results [0. However, 2D correlations are also important 
(by continuity they persist to the ferromagnetic F phase 
at saturation and to the spin liquid phase above T/y, both 
show 2D character) and give rise to an asymmetric distri- 
bution of scattering weight around [0|3). Linear 
field dependence of the incommensuration of the longitu- 
dinal spin correlations is a signature of exclusion statistics 
for the spinon quasiparticles in the HAFC, and the ob- 
served similar behavior of the incommensuration in the S 
phase dominated by the ordering of the longitudinal spin 
components suggests that exclusion statistics are impor- 
tant for the quasiparticles in CS2CUCI4. The existence 
of a modulated continuum upper boundary indicating 
phase space restrictions for paired states also supports 
this conclusion. Susceptibility measurements show 
no evidence of a phase transition between the spin liquid 
behavior in zero field above TAr=0.62 K and the disor- 
dered phase found for fields greater than 2.1 T down to at 
least 35 mK.This indicates that fields applied within the 
ordering plane stabilize the fractional spin liquid state. 

In conclusion we have studied the ground-state and 
excitations of the frustrated quantum antiferromagnet 
CS2CUCI4. This material has a 2D Hamiltonian interpo- 
lating between the square, triangular and ID Heisenberg 



antiferromagnets, and shows: (1) very strong quantum 
renormalizations indicating the importance of fluctua- 
tions; (2) continua in the dynamical correlations demon- 
strating fractional excitations; (3) very large field-driven 
incommensuration and disorder effects from inplane fields 
showing that exclusion statistics are important; and (4) 
stabilization of a spin liquid ground state by inplane 
fields. We believe new theoretical work is needed to ex- 
plain these findings. 

Full details of the analysis and extensive results from 
other related experiments on CS2CUCI4 will be given in a 
forthcoming publication |l^]. We would like to thank 
M. Eskildsen, M.A. Adams and M.J. Bull for techni- 
cal support and we acknowledge very useful discussions 
with D.F. McMorrow, R.A. Cowley, F.H.L. Essler, A.O. 
Gogolin and M. Kenzelmann. ORNL is managed for the 
US DOE by UT-Battellc, LLC, under contract DE-AC05- 
00OR22725. A.M.T. is grateful to Isaac Newton Institute 
and Trinity College Cambridge for kind hospitality. 
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FIG. 1. (a) 2D couplings in CS2CUCI4: Strong bonds J 
(heavy lines) and smaller frustrating zig-zag bonds J' (thin 
lines) . (b) and (c) Magnetic phase diagram in a field along the 
c and a axes, respectively. Symbols show the boundaries of the 
various phases described in the text measured using neutron 
scattering (squares) |llUl^ and susceptibility (circles) 
Solid curves are a guide to the eye, and dashed line indicates 
crossover to paramagnetic behavior, (d) Spin rotation in the 
(6, c) plane. Heavy arrows are spin vectors and the circle indi- 
cates the spin rotation upon translation along the fe-axis. (e) 
Incommensuration e vs. field along c (f) e vs. field along 
a (solid line is a guide to the eye, solid circle is from ||l3|). 



FIG. 2. (a) Dispersion of the magnetic excitations along 
the 6*-axis {T—0.1 K, zero field). Filled symbols (triangles 
from 1^) show the main peak in the measured lineshape 
and solid line is a fit to the principal spin- wave mode ^{k) 
of Hamiltonian (^) (dashed and dash-dotted lines show the 
corresponding dispersion of the other two secondary modes, 
w"*" and u)~ respectively, as described in the text). Typical 
scan trajectories are shown by the light shaded regions la- 
belled A-D (the line thickness is the wavevector averaging) 
and the measured data are shown in Fig. ^A-D. Open cir- 
cles show the experimentally estimated upper boundary of 
the continuum and upper (heavy) dashed line is a guide to 
the eye. Dotted area indicates extent of magnetic scattering, 
(b) Intensity measured along scan D above, in the cycloidal 
phase at 0.1 K (open circles) and in the paramagnetic phase 
at 15 K (solid circles). Data points are the raw counts and 
dashed line shows the estimated non-magnetic background. 
Solid lines are guides to the eye and the horizontal bar indi- 
cates the instrumental energy resolution. 



FIG. 3. Magnetic inelastic scattering measured at 0.1 K 
along the light shaded regions labelled A-D in Fig. ^(a) . Top 
axis shows wavevector change along scan direction. Count- 
ing times were typically 35 hours at an average proton cur- 
rent of nOfiA. Solid lines are fits to a modified two-spinon 
cross-section (see text). Vertical arrows indicate estimated 
upper boundary. Dashed lines show predicted lineshape 
for polarised cycloidal spin waves and the dark shaded re- 
gion (shown only in D for brevity) indicates the estimated 
two-magnon scattering continuum. All calculations include 
the isotropic magnetic form factor of Cu"*"^ ions and the con- 
volution with the spectrometer resolution function. 
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